Density functional theory for freezing transition of vortex-line liquid with periodic 

layer pinning 
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By the density functional theory for crystaUization, it is shown that for vortex lines in an un- 
derlying layered structure a smectic phase with period m = 2 can be stabilized by strong layer 
pinning. The freezing of vortex liquid is then two-step, a second-order liquid-smectic transition and 
a first-order smectic-lattice transition. DFT also indicates that a direct, first-order liquid-lattice 
transition preempts the smectic order with m, > 3 irrespectively of the pinning strength. Possible 
H — T phase diagrams are mapped out. Implications of the DFT results to the interlayer Josephson 
vortex system in high- Tc cuprates are given. 

PACS numbers: 74.25.Dw, 74.25. Qt 



It is now well established that the Abrikosov vortex lat- 
tice in type II superconductors 0] melts via a thermody- 
namic first-order transition J2||. However, much less con- 
sensus has been reached on the possible phases and melt- 
ing process of interlayer Josephson vortex lattice where 
layer pinning is essential. A Lorentz-force independent 
dissipation was found in Bi2Sr2CaCu208-i-y at H = bT 
parallel to the a6-plane [l] , which ispossibly a signature 
of a Kostlitz-Thouless (KT) phase [lH0. On the other 
hand, transport experiments conducted in YBa2Cu307_5 
suggested a second-order softening of the vortex lattice 
Q, which can be accounted for by a smectic phase [g. 
A recent computer simulation found continuous meltings 
and a KT-type intermediate phase above a multicritical 
field ^Jflj. 

As a general theory for crystallization, the density 
functional approach was formulated by Ramakrishnan 
and Yussouff lllj] (See also 113] )■ The basic idea of the 
density functional theory (DFT) is to describe the lattice 
phase and the freezing transition by the direct pair corre- 
lation function (DPCF) in liquid phase, which in turn can 
be evaluated from macroscopic quantities such as the av- 
erage density and temperature. At the clean limit, both 
the first-order nature of the vortex-line liquid freezing 
and quantitative aspects such as the Lindemann num- 
ber have been successfully captured by DFT [l^. Later 
on, DFT was used to investigate vortex-line systems with 
point-like |l4| and columnar 15] pinning centers. 

In the present work, we apply DFT to explore the freez- 
ing of vortex-line liquid subject to periodic layer pinning 
[H [13. As it is hard to evaluate the DPCF with fi- 
nite interlayer Josephson coupling f^, which is essen- 
tial in the present case, we will not try to draw the de- 
tailed phase diagram. Instead, we adopt the DPCF and 
the layer pinning strength as parameters to explore the 
possible freezing processes. This approach turns out to 
be fruitful: A vortex smectic phase with more vortex 
lines in every other layers (period to = 2) is observed at 
strong layer pinning. When the smecitc phase is present, 
the freezing of vortex liquid is two-step, a second-order 



liquid-smectic and a first-order smectic-lattice transition. 
A direct, first-order liquid-lattice transition preempts the 
smectic order with ttt, > 3 irrespectively of the pinning 
strength. Possible topologies oi H — T phase diagrams 
are figured out. 

The free energy of a vortex-line system measured from 
the uniform liquid is given by (12;, J_3, J_8j 
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where Sp(r) = p[v) — po with po the average areal vortex 
density and C(r — r',n— n') the DPCF. The summations 
are taken over vortex segments along the magnetic field 
(-B]|y), and the integrals over the transverse directions x 
and z, with zjjc axis. Vp is the pinning energy per vortex 
segment. 

By means of the variational calculus, the free-energy 
density per vortex segment of a vortex lattice is found as 
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with Su = l/po the area of real-space unit cell (u.c), 
Nxz the number of unit cells and A^j^ the number of vor- 
tex segments, the summation over the reciprocal lattice 
vectors (RLVs), p{r) = pg + ^j^pKexp(iK • r), and 
C(k, 0) = /d^r^„ C(r, n)exp(— ik • r). The chemical 
potential A is given by 



■Su Ju.c. 



with the integral over the unit cell. The lattice phase is 
characterized by the non-vanishing Fourier components 
of the vortex density, which thus serve as the order pa- 
rameters for the freezing transition. Following previous 
studies |il2. 13] , wc adopt a few order parameters in the 
present work confirming that more will not change the 
conclusions qualitatively. 

Because of the linear coupling with the vortex density, 
the layer pinning Vp appears only in the chemical poten- 
tial A in the free-energy density. The layer pinning is 
effective when the pinning wavenumber q coincides with 
one of the RLVs. As the DPCF C(K, 0) decreases quickly 
as going to higher shells of RLVs, the layer pinning is 
most important when its wave vector coincides with one 
of the primitive RLVs. 

For weak magnetic fields, vortices reside in a periodic 
subset of block layers with period ms where s is the layer 
spacing. For large m, the pinning energy term in Eq.(3) 
decouples from the lattice-order term and only modifies 
the measure of the integral of the lattice potential. In 
such a case the layer pinning is irrelevant to the phase 
transition. The freezing process is then essentially the 
same with the Abrikosov lattice with the scaling by the 
anisotropy parameter ll.9| . It takes numerics shown be- 
low to see to what extension that the large-7n limit gov- 
erns the freezing phenomenon . 

Vortex segment and layer pinning energy- When a single 
vortex line is placed in a periodic layer pinning potential, 
a kink (kernel) presumes a metastable state in which a 
part of the vortex line cross the pinning barrier and reside 
in the neighbor energy valley. The crossing takes place 
within the length scale L„aii — s-\/7eo/4C/p, where 7 = 
Ac/Aab, eo = (<?^o/47rAab)^, Up the pinning ene rgy per unit 
length |l,[i3|- The pinning energy is Vp ~ Sy^'yeoUp, and 
the kink length defined as the separation between two 
crossings is Lkink = iwaiiexp(2/3Vp) [3]. The physical 
choice of the vortex segment in the free-energy functional 
(1) is therefore the kink with the pinning energy per kink 
Vp. The above estimates are physically meaningful when 
ikink is sufficiently larger than Lwaii, or equivalently f3Vp 
is large enough. Numerically, for YBa2Cu307_5 with 
Tc = 92K, K = 100, 7 = 8, Aa6(0) = lOOOl and s = 
I2A, one has Lwaii — 10s, £kink — 60s, and pVp ~ 0.8 
at T = 88K. As temperature decreases, (3Vp increases 
linearly to pVp ~ 4.3 at T = 80K, Lwaii decreases slowly 
and Lkink increases exponentially. For Bi2Sr2CaCu208+y 
with Tc = 90K, K = 100, 7 = 150, A„fc(0) = 2000l and 
s = I5A, one has /3Vp ~ 0.2 at T = SAK, which increases 
linearly to pVp ~ 0.3 at T = 78K. 

Lattice unit cell- The vortex lattice structure in the pres- 
ence of layer pinning is nontrivial. In an anisotropic 
material without layer modulation, there are two sta- 
ble directions of triangle vortex lattice, which are de- 
scribed by the following two unit cells J3: Type A: ai — 

CV7x,a2 - ^x+ f^z; Type B: ai = ^z,a2 = 
^^x + 2^z, with C^ = 2(j)o/V3B. The associated 
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FIG. 1: Real-space unit cells commensurate with layer mod- 
ulation. The X direction is rescaled with the anisotropy pa- 
rameter 7. A and B are for the type of unit cell, and the 
number denotes m defined in text. 
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In DFT the reduction of free energy of crystallization 
is evaluated by the DPCF of liquid, which is isotropic 
after the space-rescaling with the anisotropy parameter. 
Just above the freezing, there appears a sharp peak in 
the DPCF at fco = 47r/\/3C associated with the uniform 
vortex density of the rescaled, isotropic liquid. The re- 
duction of free energy is maximal when the RLVs of a 
candidate lattice coincide with fco in magnitude. It is 
easy to see that the two lattice structures listed above 
match this condition. 

On the other hand, the vortex lattice is requested to 
be commensurate with the underlying layer modulation 
in DFT, which is reasonable for layered superconduc- 
tors like high- Tc cuprates. This makes the pinning en- 
ergy equivalent for different lattice structures. There- 
fore, at the following sequence of the magnetic fields 
Hi = y3</)o/27s2, H2 = </.o/2y37s', ^3 = V^hl^is'^. 
Hi — (/)o/ 6-^/375^, and H5 = 0o/8a/37s^, ..., the vortex 
lattices depicted in Fig. 1 resume minimal free energy. 
The stable lattice structure at magnetic fields nearby 
those listed above should be the same with those in Fig. 
1 except for that the vortex separation in the x direc- 
tion is adjusted in order to accommodate the vortex den- 
sity. The ground-state lattice structure changes when 
two freezing DPCFs coincide. The lattice structure with 
a smaller freezing DPCF is achieved at a higher temper- 
ature and is stable for the given magnetic field. Possi- 
ble first-order transition between different lattices upon 
sweeping the magnetic field at fixed temperature was ad- 
dressed in Ref.|20l|. 

Freezing transitions- We now proceed to investigate the 
freezing process and nature of phase transition when the 
lattice structure is given and the DPCF is swept. This 
corresponds to reducing temperature at given magnetic 
fields. We will focus on the sequence of magnetic fields 
listed above, where the vortex lattice are naturally com- 
mensurate with the layer structure. Although analyses 
on other magnetic fields are technically cumbersome as 
C(K, 0)'s assume different values even on the same shells 
of RLVs, the possible phases and nature of phase transi- 
tion should be the same. As read from Eqs. (2) and (3), 
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FIG. 2: First-order freezing into the lattice structure A2 
for /?Vp = 0.5, with C'^^ = 0.3C'^>. Thermal fluctuations 
(it^) and (u^) are normalized by (7s)^ and s^ respectively. 
\ul) = oo for poC'^) < 0.7382. 
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FIG. 4: poC'^^-' —f3Vp phase diagram for the lattice A2. Inset: 
mean squared displacements, normalized by (75)'^ and s'^ in 
the X and z directions respectively, evaluated along the first- 
order phase boundary. 
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FIG. 3: Two-step freezing, a second-order and a first-order 
transitions, into the lattice structure A2 and an intermediate 
smectic phase for /3Vp = 1, with C*^' = 0.3C'^'. Thermal 
fluctuations (u^) and {u^} are normalized by {-ys)^ and s'^ 



respectively. 



CX3 for poC'^> < 0.6856. 



freezing transitions for different anisotropy parameter 7 
are governed by the same set of C(K,0)'s as functions 
of Vp and m. This corresponds to the scaling property 
addressed in Ref . [l9| . apart from the layer pinning effect. 

For Al in Fig. 1, the present system is equivalent 
to colloids under laser radiation addressed in Ref. [l8l| . 
As the pinning energy increases, the freezing switches 
from first order to continuous at a tricritical point at 
PVp ~ 0.212 and poC^^) ~ 0.748. The notation with C^*) 
referring to C(K,0) for K on the i-th shell of RLVs is 
adopted in the present Letter. 

For A2, we use four order parameters pK where Kl = 
b2, K2 = bi, K3 = 2b2, and K4 = 2bi with ba = tt/sz. 
A finite pK3 characterizes the modulated liquid. The 
freezing is first order at low pinning energy such as 
f3Vp = 0.5 depicted in Fig. 2, where all order param- 
eters except for the trivial one set up discontinuously. At 
large pinning energy, a phase characterized by pKi > 



and pk2 = Pk4 = appears via a second-order phase 
transition as shown explicitly in Fig. 3 for pVp = 1 . The 
vortex density modulation is associated with a wave num- 
ber half of the underlying layer structure, corresponding 
to more vortices in every other layers; the system be- 
haves as liquid in in- plane directions. This is the smectic 
order addressed in Ref.^], and the liquid-smectic tran- 
sition is in the Ising universality class since m — 2. As 
DPCF increases, the full lattice order appears via a first- 
order transition. The poC^^' — (3Vp phase diagram for 
A2 is presented in Fig. 4, with a multicritical point at 
poC(i) ~ 0.714 and (3Vp ~ 0.712. This is one of the 
central observations of the present DFT study. 

Similar analysis has been performed for A-type lattice 
withm > 3, using C(2) = 0.5C<^^\ C^^) = 0.3C<-^\ C^^) == 
O.IC^^^ and so on. In contrast with Al and A2, we cannot 
find the smectic phase and continuous transition up to 
huge pinning energy accessible in numerics. This DFT 
prediction is different from the analysis in Ref. 8] . 

For unit cell of type B, the crystallization is always one 
step and first order. In this case no primitive Bragg peak 
matches with the layer modulation, and thus the layer 
pinning is much ineffective compared with the lattices Al 
and A2. As the result, the six primitive Bragg peaks ap- 
pear simultaneously which is associated with a first-order 
transition. The first-order melting transition observed in 
Monte Carlo simulations for B = (/)o/32s^ and 7 = 8 
[3, U^ is in accordance with the DFT, since the lattice 
structure is Bl. In principle, there might be a smectic 
phase associated with the two Bragg spots K = (0, ±tt/s) 
on the second shell of RLVs for B2. However, numerical 
analysis reveals that this situation is not realized since 
the layer pinning is at the 6th shell of RLVs and the 
interaction with those on the 2nd shell is very weak. 

Phase diagram- Although the detailed freezing curve 
cannot be drawn without calculations on DPCF at given 
magnetic fields and temperatures, we are able to catego- 
rize possible H — T phase diagrams based on the above 
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FIG. 5: Possible H — T phase diagrams for vortex lines in 
periodic layer pinning potential from DFT, with L, X and 
S abbreviating liquid, lattice and smectic respectively. The 
thick (thin) curves denote first(second)-order phase bound- 
aries. Multicritical points are marked by open circles. The 
horizontal phase boundaries are drown only for B1-A2 (A1-A2 
in (d)) and A2-A3 transitions. 



analyses. Figure 5(a) corresponds to a system of weak 
layer pinning, in which freezing is always first order. 
The two horizontal phase boundaries are between lat- 
tices structures B1-A2 and A2-A3 ,20J. Note that simi- 
lar phase boundaries may be observed at lower magnetic 
fields, and that the lattices Al and Bl can be tuned into 
each other through crossover. In Fig. 5(b), the freezing 
transition into vortex lattice Al is of second order at high 
magnetic fields where the freezing temperature is low and 
the layer pinning is strong enough. At even stronger layer 
pinning, a smectic phase appears and transforms into 



lattice A2 upon cooling as shown in Fig. 5(c). As the 
magnetic field increases to the regime where the lattice 
Bl presumes the ground state, the smectic phase shrinks 
to zero. It is possible practically that the lattice Bl is 
unstable comparing with other lattice structures, result- 
ing in Fig. 5(d). No further complex phase diagram is 
expected from DFT. 

Discussions- Based on the DFT analysis and experi- 
mental phase diagrams for H < lOT [3, l21|, the m = 
2 smectic phase is expected for YBa2Cu307_5 around 
H = H3 ~ 39T. No smectic phase is plausible for 
H < H4 ~ (/)o/6v37S^ ~ 17T. Layer pinning may be 
too weak in Bi2Sr2CaCu208+y to stabilize the smectic. 

Since DFT only treats the RLVs, it is unable to de- 
scribe quasi long-range orders. A KT phase is therefore 
out of the scope where DFT can reach. With this draw- 
back of DFT in mind, one may notice that the tricriti- 
cal point for Al lattice in Fig. 5, which should exist in 
both YBa2Cu307_5 and Bi2Sr2CaCu208+i;, corresponds 
to the multicritical point above which an intermediate 
KT phase is observed in computer simulations jlOl |. 

As displayed in the inset of Fig. 4, thermal fluctua- 
tions are more anisotropic than the anisotropy parame- 
ter 7, which is clearly caused by the layer pinning. This 
anisotropy exists irrespectively of the existence of smec- 
tic phase, and thus cannot be taken as a precursor of 
a partial melting |22| |. In presence of layer pinning the 
Lindemann numbers are not constant any more. 

One may expect that the results derived above apply 
to two-dimensional (2D) systems jl2| with periodic line 
pinning. It is noticed however that thermal fluctuations 
are more important in 2D as addressed in Ref. |23| . 
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